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TECHNICAL NOTE 2959 

THEORETICAL INVESTIGATION OF THE SUPERSONIC LIFT 
AND DRAG OF THIN, SWEPTBACK WINGS WITH 
INCREASED SWEEP NEAR THE ROOT 
By Doris Cohen and Morris D. Friedman 


SUMMARY 


Formulas are derived by the use of linearized theory for the lift 
and drag due to lift, at supersonic speeds, of thin, flat wings having 
a discontinuity in the leading-edge sweep, the inhoard portion of the 
leading-edge being very highly swept and the outboard portion less so. 
Examples are presented to show the effect of the bend in the leading 
edge on the pressure distribution. The lift-curve slope and drag have 
been calculated for several families of wings, all with straight trail- 
ing edges. For two typical plan forms, the aerodynamic-center location 
has been calculated through a limited range of supersonic Mach numbers. 

The over-all characteristics of the wings studied show little 
effect of the concentration of sweep near the root, in the absence of 
thickness and viscosity, but appear to be determined primarily by the 
sweep of the outer portion. However, there is a shift of the lifting 
pressure away from the central portion of the wing and toward the lead- 
ing edge of the outer portions . 'In most cases , there will also be a 
region of high lift around the trailing edge of the root section. As 
a result, the aerodynamic center is general 1 y farther back than on com- 
parable conventional wings. In the limited calculations made, no shift 
of aerodynamic -center location with Mach number was observed. Similar 
effects are to be expected on the wing of a wing -fuselage configuration 
of similar plan form. 


INTRODUCTION 


Considerable interest has been shown in the 'use of wing plan forms 
incorporating a region of increased leading-edge sweep, similar to a 
large fillet, near the wing root. The forward extension of the root 
chord can provide both increased wing depth, with its structural 
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advantage, and a useful and accessible volume for the installation of 
fixed equipment, fuel, air scoops, or complete power plants, depending 
upon the airplane size and purpose. Such a wing plan form should he 
particularly useful for an all -wing airplane because of the greater 
flexibility allowed in the chordwise location of the major weight 
items. 

The use of a so-called "cranked" or "crescent" plan form also has 
several potential aerodynamic advantages. 1 Theoretical considerations 
indicate that the low-speed stability and control characteristics should 
be improved and the high-speed wave drag reduced by increasing the sweep 
near the root and decreasing it near the tips. Early attempts to real- 
ize these advantages resulted in the experimental wings reported in 
references 1 and 2. The wings did, in fact, show good maximum lift and 
aileron effectiveness at low speeds, although no extensive comparison 
was made with conventional swept wings . Reduction in form drag at high 
subsonic speeds was achieved because the wing chord was increased with- 
out increasing the thickness. 

As far as is known at this time, no investigation has been made 
into the supersonic characteristics of the cranked wing. The present 
report considers the theoretical forces and moments on such wings due 
to angle of attack, for speeds greater than the speed of sound. In 
order to reduce the number of variables, the wing plan form has been 
kept as simple as possible. In every case, the trailing edge is unswept. 
A single bend is introduced Into the leading edge at various locations, 
and various amounts of sveepback are specified for the inner and outer 
portions. The sweep of the inner portion of the wings Is always high 
enough so that the forepart of the wing lies Inside the Mach cone from 
its apex. Finally, the wings considered have tips either tapered to a 
point (type A of fig. l) or cut off in the streamwise direction (type B) . 
The distribution of lift, the total lift, and the drag due to lift are 
derived using linearized supersonic wing theory, and some investigation 
of the aerodynamic-center location is made. 

The applicability of the results to the all -wing airplane with 
thickened and lengthened root section has been mentioned. The results 
should offer, also, at least a qualitative indication of the distribu- 
tion of pressures to be expected on the wing of a wing -fuselage com- 
bination of similar plan form, particularly if the fuselage has a 
flattened cross section. For information bearing on this point, the 
reader is referred to a paper by Morikawa (ref. 3 ) ■ 


1 Both structural and aerodynamic considerations behind the selection of 
the crescent wing are discussed in an article by R. S. Stafford in 
the British publication, The Aeroplane, Jan. 2, 1953> PP* 6-7. 
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METHOD OF ANALYSIS 
Approach to the Problem 


The method of analysis will he described in some detail, since it 
is applicable to a more general class of wing than is discussed here. 
Basically, the method consists of the superposition of solutions of the 
linearized supersonic-flow equations. For the most part, conical-flow 
solutions will he employed, in the manner suggested by Busemann (ref. 4) 
and Lagerstrom (ref. 5) • The method of Mirels (ref. 6 ) will also be 
used, however. 

Consider first wing A. Nothing that occurs behind its trailing 
edge can affect the flow over the wing itself (since small perturbations 
are assumed). Hence , it is permissible to ignore temporarily the posi- 
tion of the trailing edge, that is, to assume that the wing extends to 
infinity both laterally and to the rear . Now consider the infinite 
triangle formed by the downstream extension of the inner portions of 
the leading edge. (See sketch (a).) 
triangle flying as an isolated wing 
is, of course, known (refs. 7 > 8 , 
and 9 ). However, the accompany- 
ing upwash on either side of the 
triangle violates the require- 
ment of tangential flow over the 
outer panels of the complete wing. 

We therefore superimpose addi- 
tional solutions of the supersonic 
flow equations which have the prop- 
erty of "canceling" the upwash 
within the area of the horizontal 
plane occupied by the outer panels, 
without introducing any new verti- 
cal velocities through the region 
occupied by the center triangle. Sketch (a) 

This step determines that the flow will conform to the condition of a 
flat plate outboard of the triangle; it is further necessary to intro- 
duce a constant downwash in these regions to correspond to the angle 
of attack of the wing as a whole. This last must also be achieved 
without the introduction of additional downwash inboard of the panels . 
The net result is a change in lift on the portion of the original trian- 
gular wing within the Mach cone from the apex of each panel, which, 
together with the lift on the panels themselves, is to be calculated. 


The complete solution for such a 
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Notation 


The symbols used in this report are defined in Appendix A. In 
addition to the Cartesian coordinate system (fig. 2) with origin at 
the apex of the wing, x axis positive downstream and z axis positive 
upward, a conical coordinate system in which the location of the origin 
is arbitrary will be used. Observed in the xy plane, the character- 
istic variable of such a system is the inclination of the line from the 
chosen origin x±,y± through the variable point x,y, the reference 
being the free -stream direction. For supersonic wing theory, the ratio 

of this quantity to the tangent of the Mach angle 1 /Vm 2 - 1 is most 

significant. The conical variable 
(see sketch (b)) will therefore be 
chosen as 

y - y± , — 5 

t = 0 (3 = v M 2 - 1 

X - Xj_ 

with subscripts to indicate the 
location of the origin. . For brev- 
ity, a line making the angle 

tan -1 t/p with the stream will 
frequently be referred to as "the 
ray t." The lines t = ±1 are 
Mach lines, regardless of the 
origin. 

In keeping with the foregoing coordinate system, the sweepback of 
the leading edge of the wing is expressed in terms of similar parameters: 

mi= (3 times the tangent - of the semiapex angle of the wing, or the cotan- 
gent of the angle of sweep of the inner portion of the leading edge 

m 2 = (3 times the cotangent of the angle of sweep of the outer portion 

The linear dimensions of the wings are defined in figure 2 and 
Appendix A. 



Boundary Conditions for Elementary Solutions 


The upwash field of the inside triangle can be written from refer- 
ence 8, in terms of the conical variable t Q = (3 y/x (see fig. 3 and- 
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Appendix A) as follows: 


v A^o) = Sr 


1-t 


t 0 2 -m 3. 2 


2 E(<P, a/ I-K! 2 ) 


s mi < t Q < 1 


(1) 


where E is the elliptic integral of the second kind, with the argument 



and u has the value 


miVa 

Uo = 

PECVl-ni! 2 ) 


(2) 


(3) 


The solution required to cancel this field is constructed by the 
superposition of a family of solutions with parameter t 0 , each solution 
having the following properties: 

1- In the infinite sector between any ray By = t 0 x and the 
outer segment of the leading edge (hatched portion, fig . * 0 , the 
vertical velocity is constant. 

2. Between the ray gy = tox and the Mach line crossing the 
wing from the intersection xi,yi of the ray with the leadiDg edge, 
the vertical velocity is zero, (it is assumed for the moment that 
the Mach line does not cross the leading edge of the opposite panel.) 

3- Along the Mach line just described, all the velocity compo- 
nents of the cancellation field vanish. 

2. Ahead of the leading edge (on the right, in fig. 2) the 
streamwise component of the perturbation velocity is zero . 

The foregoing conditions define a conical flow field (ref. 2), that 
is, one in which the velocity components are constant along any ray from 
the origin or apex - in this case the point xi,yi. The solution is 
therefore expressible as a function of the single variable 



x - X3. 


the ratio of the slope of a ray of the field to the inclination of the 
Mach lines. 
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The conditions specified are the same as for the flow field of a 
deflected, triangular leading-edge flap, and the required solutions haye 
been given in that connection in reference 5* In the sections to 
follow, these solutions will be used to build up the required wing plan 
forms by the method of superposition, and the resulting lift and drag 
will be calculated. Because of differences in detail between the calcu- 
lation of the lifting pressure at any point and of the integrated lift 
and drag, the remainder of the analysis will be presented with specific 
application to each of these problems. 


LIFT DISTRIBUTION 


Suppose the incremental lifting pressure coefficient associated 
with each elementary solution is 


AC p (x,y,t 0 ) = - ^ | W(t 0 )P(t 0 ,t 1 ) (5) 

in which W is the constant vertical velocity on the sector and V the 
stream velocity. Then the total increase in local lift at a point x,y 
due to canceling the upwash on one side of the wing is 


8 [ 

rt L V 


B?Oi,t m ) 



<3w A (t Q ) 

dt Q 


P(t Q ,ti) 



(6) 


where tm is the value of 


ti when to = mi, or 


tm - 3 


y 

x 


Sl 


Cl 


(7) 


and T 0 (x,y) is the value of t 0 such that 

ti(x,y,r 0 ) = -1 


Since 


xi 


m. 


mi 

Cl 


m 2 - t 0 


( 8 ) 
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and 


3yi = toxi 

setting tx (eq. (ij-)) equal to -1 gives 


(9) 


T oU>y) = 


m 2 (x+ jBy) - (m £ -m 1 )c 1 

( x+ 0y) + (m 2 -mx) ex 


( 10 ) 


Also needed for equation (6) is the derivative of 


Wa’ ( t Q ) = -mx Uq 


^ l-t Q 2 
,2 V 3 / 2 


(t 0 2 -mx 2 ) 


In addition to canceling the downwash, it is necessary to take into 
account the effect of the incidence a of the panel. This effect is 
obtained from equation (5) with W = -aV and to = so that the total 
lifting pressure induced at a point x,y by the addition of the panels 
to the inner triangle is 


AC 


P 


- 2 


w ; 


a+ • 


v (mi) 


P r o 




w. 


(to) 


V 


■S>(to,tx) dt G [ (11) 


If the point x,y is on the inside triangle, the foregoing lift 
is added to the lift associated with the original triangular -wing 
solution: 


Cp(x,y) 

A 


4m.iUoX 

V>/mi 2 x 2 -p 2 y 2 


( 12 ) 


and, if it lies within the Mach cones from the apexes of both left- and 
right-hand outer panels, ACp(x, -y) must also be added. For wing B, 
tip corrections will also be derived, so that, for a point near the 
center of the trailing edge, the complete expression for the pressure 
might be 


C p (x,y) =Cp A (x,y) + ACp(x,y)+AC p (x,-y) + (AC p ) tip (x,y)+ (ACp) (x, -y) 

(13) 
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Elementary Solutions 


At this point it "becomes necessary to specify whether the sweep of 
the Mach lines is greater or less than the sweep of the outer segment 
of the leading edge, that is, whether ma is greater or less than one. 
For mad, the solution satisfying the specified conditions is given 
"by equation (6.6) of reference from which is obtained 


[]P(to,ti)] m = m 2 PiCto^tx) - t Q P2(to,t!) 


( 1 ^) 


with 


Pi(t 0 ,tx) 


/(m 2 - t Q ) (1+ t x ) 


1+maV (1 - t Q ) (m 2 - t x ) 


( 15 ) 


and 


x y(l + t 0 )(m 2 - t a ) + J (1 +ti)(me - t 0 ) 

P 2 (to,t!) = log — — 

Vl-to 2 V(l+m2) Jt 0 - t x f 

These functions and the resultant pressure are plotted against ti in 
figure 5 for typical values of t 0 and m 2 . Note that the pressure is 
infinite along the subsonic leading edge and also when t x = t G , the 
hinge -line location in the equivalent control -surface problem. 

For m 2 >1, the solution is identical with that on a vertically 
symmetrical wedge, with the lower-surface values changed in sign. The 
flow over the wedge, which is required to have a constant vertical 
velocity in the region between the leading edge and the extension of 
the ray t Q across the wing (sketch (c)), may be constructed by sub- 
tracting from the flow field of a wedge, bounded by the leading edge 
and a line in the downstream direction from Xi,yi, the velocities over 
a wedge of equal inclination, bounded by the ray t 0 and the line 
drawn downstream from xi,yi< The two component solutions may be 
obtained .from reference 5 , or more readily, from reference 10. The 
combined solutions give a pressure proportional to 
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and 


QsC'toj'ti) - 




■ log- 


J (l+ t Q ) (1 - tl) + J { 1 + tl) (1 ~ "bo) 

^2|t 0 - ti f 


(19) 


These functions and the resultant pressure coefficient are plotted in 
figure 6. 

Formulas and Lift Distribution - Leading Edge Entirely Subsonic 


Formulas for If 3P(t 0 ,ti) (eq. (l4)) is substituted in equa- 

tion (11), the latter may be integrated (ref. ll) to give the following 
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result: 


AC p (x,y) = - ~ -j m lg (t,T 0 )G n - 


m 2 (l + t) 


2m 1 (l+ m 2 ) g (t, t d ) 


[2E(k) - k' 2 K(k) ] 


where 

t 


py 

x 


k = 


(1 -m 1 )(T Q - m x ) 
(1+ mi) (t + 


g = 


(l+t) ( t q -t) 


(l + m 1 )(T 0 + m 1 ) 

k' = J 1 - k 2 


(20) 


and 


G x = K(k) 


Z(Si,k) 


L mi - t g J t z - mi 2 


„ , /(1+ mi)(t+ mi) 

; 0i = sin / for t < -mi 

v (1 - mx) (t - mi) 


Gz '■ 


-1 


1 - mi 


1 + mi 

K(k) E(k) 

2m x 


for t = -m x 


Go = 


K(k) it A c( e s^ k ) 4 


, /(l + mi) (t + mi) 

0 g = sin / — for -mi < t < mi 

l + t 2 J J 2mi(l + t) 




-K(k) 


1 

t + mi 


z(e 4 ,k) i 

____ • 

g'J t 2 - mi 2 - 


0 4 


sin -1 / 

J 


(t- mi)(T 0 + mi) 
(t+ mi) ( t q -mi) 


for 


t >mi 


The elliptic functions E(k), K(k), Z(0,k), and A Q (0,k) are defined 
in Appendix A. 

Along the inboard Mach line y = y*(x) from the apex ci,si of the 
panel, t q reduces to mi, and AC p takes on the simple form 
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It-Uo - mi ) PsI 

V mi(l+ mg) v 2(1 + mi) (x - Ci) 


( 21 ) 


Ahead of this line the effect of the panel is of course zero, so that 
the result is a finite step in the lift distribution. It is interest- 
ing to note that, with m 2 = 0, equation (2l) gives the value of the 
decrement in lift along the corresponding Mach line from the outside 
corner formed by a streamwise tip on a simple sweptback wing (see 
sketch (d) and ref. 12). 


Along the boundary between 
the outer panel and the inside 
triangle, t equals m x , G 3 is infi- 
nite, and ACp is double -valued. 
However, both the limit of Cp +ACp 

as t approaches mi from below, 
and the limit of ACp as t 
approaches m x from above (where 
Cp =0) equal 



Cp 



8uq m i f T 0 + mi 

ftV Vt 0 2 -mi 2 i ^mi 


m 2 (l + mi) 
m x (l + m 2 ) 


(2E -k' 2 



( 22 ) 


with the modulus of the elliptic integrals as defined for equation (20) 
It will be seen from the numerical results that the lift-distribution 
curve goes smoothly through this value, although the initial lift dis- 


tribution Cp went to infinity along the line t = mi. 


If the Inboard Mach line 
from ci,si intersects the oppo- 
site leading edge as in sketch (e) , 
fictitious pressure differences 
will have been introduced, by the 
use of the conical fields, into the 
region ahead of the wing shown 
shaded in the sketch. Cancella- 
tion of these pressure differences 
is necessary in order to satisfy 



Sketch (e) 
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the condition of continuity of pressure in the stream. Mathematically, 
this is a somewhat formidable task. As the pressures involved are 
actually not very great, it will usually be sufficient to perform only 
an approximate cancellation, as follows: 

Assume ACp to be constant and equal to ACp (x*,y*) throughout 
the shaded region. Then it may be canceled, and the boundary conditions 
on the wing preserved, by a single conical field with its origin at 
x*,y*. In terms of the conical variable 


t* = 0 


y- y* 


X - X* 


the boundary conditions on the solution are simply. 


— = AC p (x*,y*); -1 < t* ^ -m 2 
w = 0; -m 2 < t* < + 1 


The required solution is obtained by a slight modification of equa- 
tion (5.15'ti) of reference 5: 


(ACp) L.E. (x ' y) = ' « AC P (x *’ y *) 


2(m 2 + t*)-(l -m 2 )(l- t*) 
cos -1 

(l + m 2 ) (1+ t*) 


2m 2 V 2(1-1112) /l - t* 


1+ m 2 


m 2 + t* J 


( 23 ) 


This term is to be added into equation (13) . 

Tip effect, wing B .- If the trailing edge of the wing is unswept, 
or is swept (either forward or back) but is everywhere supersonic, and 
if the tips are tapered to a point, the entire lift distribution due to 
angle of attack is given by the foregoing formulas. If the wing is not 
fully tapered, but has tips formed by any entirely subsonic curve, the 
lift on the portions of the wing included in the foremost Mach cones 
from the tip will be modified in a manner best calculated by the method 
of reference 6. If the tips are simply cut off parallel to the stream, 
the pertinent formula (17a) in reference 6 can be integrated after the 
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insertion of ACp from equation (20) 2 to give 


x . y -SmxUo / (l+t)(T 1 -t) r /vST , . 

(^p) ttp ( x ^y) a --7 \ g 4 ( t i) /— g 4 ( t o) /-S— 

tip rtV V (l+mi)(T 0 +mi) L -/ Ti+m! v T X -t 


+ 


m 2 (l+t) (t+mi) 


/T 0 -t 


2m 1 2 (l+m 2 )g 2 (t,T 0 ) V T !-t L 


2E(k) - k’ 2 K(k) 


(24) 


in which G 4 ( t i) is obtained by replacing t q by 


Ti = 


Ps 

x - p(s-y) 


(25) 


everywhere 
tion (20) . 
reduces to 


in G 4 , and k has the same value otherwise as in equa- 

Along the Mach line from the tip, T i equals t d and (ACp) 

* tip 


(ACp) 


tip 




km 2 UQ(l+t) (t+mi) 
m 1 itV(l+m 2 )g(t,T 0 ) 


2E(k)-k ,2 K(k) 


( 26 ) 


At the tip, the right-hand side of equation (2k) becomes indeterminate, 
but reduces in the limit to -ACp(x, s) as given by equation (20) . 


Numerical results . - The lift distribution at M = *[~2. was calcu- 
lated for a wing with semiapex angle equal to tan -1 0.375 (that is, 
almost 70° sweep of the leading edge at the root) and kS° sweep of the 
leading edge of the outer panels. With p = 1, the sweep parameters 
were therefore mi = 0-375 and m 2 = 0.9* The linear dimensions of 
the wing and the Mach line pattern are shown in figure 7* 


2 It is assumed in this investigation that the region which has been 
referred to as the "inside triangle" is in fact a triangle, that is, 
that the inner portion of the leading edge, extended rearward, inter- 
sects the trailing edge rather than the tip. Formulas for the latter 
case can be derived by including Cp in the lift to be canceled at 
the tips . A 
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The resulting lift distribution is shovm (fig. 8) in spanwise 
sections, that is, for constant values of x. Ahead of and at the 
section x = 4, the lift coefficient is merely that (eq. (12)) on the 
triangular wing formed by the inner portion of the leading edge. At 
x = 4.5, the effect of the bend in the leading edge appears at y = ±1. 
The resulting step in the lift distribution mentioned earlier is seen 
to be quite pronounced. Between the Mach lines, the lift retains the 
relatively low value associated with the slender triangle. 

At x = 5.5, the Mach lines intersect on the center line of the 
wing so that the steps do not appear individually. At x = 8, the Mach 
cones overlap and the effect of the outer panels is doubled in the 
center. The step at the Mach lines has fallen off considerably in mag- 
nitude . 

The section at x = 9 , otherwise similar to the preceding one, 
passes through the tip Mach cone. Again the discontinuity in lift 
appears at the Mach line. As in the case of the simple sweptback wing 
(ref. 12), the tip effect is so large as to destroy the major part of 
the lift within the Mach cone. 

The span loadings obtained for this wing and for the type A wing 
derived from it by cutting off the rearmost portion (x >8) are shown 
in figure 9* Also shown is the elliptical span loading of a triangular 
wing. In the case of the type B wing, the tip effect appears as a dis- 
continuity in the slope of the span load curve. The bend in the lead- 
ing edge does not, however, cause any abrupt change in the shape of the 
loading, as far as can be seen. 

Aerodynamic center .- The pressure distributions of figure 8 were 
used to calculate the aerodynamic -center locations for the same wings 
for which the span loadings were calculated. The aerodynamic -center 
location was also determined at a Mach number of 1.2 and, using slender- 
wing theory (ref. 13) , at M = 1.0. This last result Is the same for 
both plan forms since, according to slender-wing theory, the non- 
expanding part of the wing is also nonlifting. 

The lengthwise distributions of load, from which the centers of 
pressure were calculated, are shown in figure 10. The same curves apply 
to both wings up to x = 8, where the type A wing ends. 

The aerodynamic -center locations are shown in figure 11 . The aero- 
dynamic center of the type A wing remains essentially fixed at x = 6.04 
or 0.755 c Q for all three Mach numbers and thus lies considerably 
behind the center of area of the wing, which is at x = 0.710 c Q . The 
aerodynamic center of the blunt -tipped wing shifts from x = 6.04 or 
0.6o4 c Q at M = 1 to 0.682 c 0 for M = 1.2 and 1.4, but still lies 
ahead of the center of area at 0.696 c Q . These results may be compared 
with the triangular wing, for which the aerodynamic center, according 
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to linear theory, coincides with the center of area at 0.667 c Q through 
out the supersonic speed range. 


Lift Distribution - Leading Edge of Outer Panel Supersonic 


Formulas for ACp.- The substitution of the elementary solution 

(eq. (17)) for m 2 >l in the integrals for ACp (eq. (ll)) results in 
hyperelliptic integrals which cannot he evaluated in terms of tabulated 
functions, except for the special cases y = y*(x) (the inboard Mach 
line from Ci,Si) and points along the leading edge. Along the inboard 
Mach line, the result is 


AC p (x,y*) 



(1 - mi) (m 2 - mi) ftsi 
(1+ m!)(l + m 2 ) (x - ex) 


(27) 


Along the supersonic leading edge, the pressure is finite, having the 
value 


( c p) 




L .E . */ m 2 2 - 1 


w 


a + 




( 28 ) 


This latter is the two-dimensional value of the pressure coefficient, 
with the angle of attack increased by the local angle of upwash of the 
streamlines due to the projecting triangular region at the root. 

Along the outer Mach line from the bend in the leading edge, the 
pressure is theoretically infinite, since the pressure for the elemen- 
tary solution (eq. (17)) does not vanish along ti = 1 for the super- 
sonic leading edge, and the magnitude W is infinite at t 0 = mi. The 
resulting infinity is of power l/2, as was previously found by Morikawa 
(ref. 3 ) for the special case of a concave 90° corner. 

The pressure on the wing elsewhere than along the lines just dis- 
cussed must be found by numerical integration of equation (ll) . Since, 
with m 2 > 1, 3P(t 0 ,ti) is constant when t x is greater than 1 (fig. 6) 
equation (ll) may be written, for points ahead of the Mach line from the 
bend in the leading edge. 


AC p (x,y) = 


4m 2 
*/ m 2 2 -l |_ 


a + 


v a ( t 2 ) 


8 

+ * 


T o w A * (t ) 


V 


P(t Q ,t 1 ) dt Q (29) 
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in which 


(m 2 -m 1 )c 1 - m 2 (x-£y) 

t 2 = 

(m 2 -m 1 )c 1 -(x-py) 


(30) 


is the value of t Q such that ti(x,y,T 2 ) = 1. 

If the point x,y is behind the Mach line, bo that t 2 is less 
than mi, the presence of the singularity in w^' at t Q = m 1 makes 
it advisable first to integrate equation (ll) by parts, obtaining 


AC p (x,y) = - 

e it 


a P (mi, tm) + 


PC®!) W^(t 0 ) 


dP(t Q ,ti) 


(31) 


Tip effect .- If the Mach number, or aspect ratio, or both, are so 
large that the Mach lines from the apex of the wing intersect the trail- 
ing edge, as in sketch (f), the pressure coefficient in the shaded 
region is constant at the two-dimensional value 



Cp - 


4m 2 a 
J m 2 2 -l 


(32) 


Then the tip effect can be found 
as the result of the superposi- 
tion of a single conical flow 
field to cancel that pressure in 
the region indicated outboard of 
the wing tip in the sketch. The 
induced decrement in lift will 
be, from reference 5> 


(AC p ) (x,y) 
e tip 


-4m 2 a 

— - cos-i 

1 ttj m 2 2 -l 


m 2 + t B + 2m 2 t s 

t S - 


(33) 


with 


"fcs 


= P 


y - s 

X - c 2 


( 3*0 


If a part of the tip section lies within the upwash field of the 
forward part of the wing, the effect on the lift distribution may again 
be calculated by Mirels' method. Inserting equations (29) and (31) in 
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equation (17a) of reference 6 and inverting the order of integration 
yield the following expressions for the tip-induced decrement to the 
pressure coefficient: 

For points in region ABDE (see sketch (g)) 



with 

i = (x-c 2 ) - p(s-y) 

(x-c 2 ) + p(s-y) 
m 2 - 1 
m 2 + 1 

(m 2 -l) ps - m 2 i 
(m 2 -l) c 2 - § 

H(t 0 ) = tan -1 





Sketch (g) 




k(t 0 ) 


(m 2 -l) 


Ps - c 2 t Q 
m 2 - t Q 


CD ^ 
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For points in region BDC, 


( ACp) tlp (x,y) - - ; 


m 2 


4 


mo -1 


o nH(li) 

a HCmi) + - J w A (t 0 ) dH(t 0 ) 


+ 


m i u o / (l+t)(r 1 -t) 1 

v 4 (1+mi) (T 1 +m 1 ) J 


(35b) 


(See eqs . (20) and (2^+) . ) 


Numerical re stilts .- The lift distribution has been calculated for 
two spanwise sections of the previously considered wing, with 70° and 
of sweep, at a Mach number of I. 667 , so that = 0.5 and m 2 = 1.2. 
The Mach line pattern and the resulting lift distributions at x = 5 
and x = 10 are presented in figure 12. 


LIFT AND DRAG 

Method of Calculating Total Lift 


In applying the conical-flows method, the simplest procedure for 
calculating the lift is to determine the lift due to each conical field 
separately by integrating the expression for the pressure with respect 
to the conical variable of that field. Thus the lift on the inside 
triangle Is found first, neglecting interference effects from the panels. 
The result may be written 



miu 0 
2n 

pv 


c 0 


(36) 


(it is assumed, as before, that the leading edge of the triangle does 
not intersect the tip of the wing.) 

An expression for the incremental lift due to each elementary can- 
cellation field is derived by integrating equation (5) with respect to ti 
over the proper area, which varies with the constant of the field, t Q . 

The results are then combined in precisely the same way, except for the 
limits of integration, as in obtaining the expression for Cp. Thus, if 
the integrated lift per elementary cancellation field is given by 
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AL(t 0 ) 

1 


-4 w(t Q ) 

« (3V 


IL (t Q ) 


the total lift on and due to the outer panels is given by 


(37) 


AL _ 8 

q. up 


w. 


a + 


A 


(mi) 


V 


]L (mi) + 


T 0 w A > (t Q ) 


V 


IL(t 0 ) dt c 


mi 


(38) 


which corresponds to equation (ll) for the incremental pressure. (A 
factor of 2 has been introduced to provide for both right- and left- 
hand panels.) The upper limit 


To = 



(39) 


of the integral is the value of t G at the tip of the leading edge. 
If the tip is outside the Mach cone from the origin, the integral is 
of course taken only to t Q = 1. 


In the case of a wing of type B, it is necessary to take into 
account the presence of the tip, not only as it affects the area of 
integration, but also as it reduces the pressure at points within the 
tip Mach cone. The latter effect, previously calculated by Mirels' 
method, is more readily computed by the conical -flows method when the 
integrated lift is required. 

If it is assumed as before that the line t Q = m x does not inter- 
sect the tip, it _Js necessary to consider the tip-interference effect 
in the derivation of IL(t 0 ) only. For each elementary upwash- 
cancellation field with apex at x x ,y x (fig. 13 ) , the lifting pressure 
varies conically according to equation (l4) or (17) between the Mach 
line t x = -1 and the leading 
edge t x = m 2 , extended to infin- 
ity. The lift outboard of the 
tip must, however, be brought to 
zero. For a single conical field, 
this is readily accomplished by 
the superposition of other coni- 
cal fields in a manner entirely 
analogous to the cancellation of 
upwash. The necessary formulas 
and the details of the procedure 
may be found in reference 12. 

If the Mach line t x = -1 
intersects the leading edge or 
the tip of the left-hand wing 
panel (as in sketch (h)), 
rather than the trailing edge. 



Sketch (h) 



20 


MCA TN 2959 


similar corrections to IL(t 0 ) should be made for the lift canceled at 
those boundaries . These effects are, however, mathematically very 
complex and in the calculations to be reported were incorporated, in 
an approximate way only, at the end of the computations. Cases in 
which these corrections might be any considerable fraction of the total 
lift were excluded from consideration. 


Men the integration of AC p (t 0 ,ti) is performed over the area 
ABCDEF of sketch (h), it is found to result in terms identical in value 
but opposite in sign to terms in the previously derived tip correction. 
It is most convenient, therefore, to use for IL(to) ia equation ( 38 ) 
the value from the combined or net lift. The resulting formulas are 
listed for the various cases in Appendix B. 


Unfortunately, in the general case, substitution of IL(t 0 ) in 
equation (38) results in incomplete elliptic integrals of the third 
kind, which are as yet untabulated and exceedingly tedious to compute . 
There is nothing to be gained, therefore, by performing the integration 


analytically, and the equation is 



left for numerical evaluation. 

In the special case shown in 
sketch (i), in which the Mach number 
is so high that the entire tip is 
outside the Mach lines from the apex 
of the wing, the elliptic integrals 
become complete, and the resulting 
formula is readily computed with the 
aid of the tables of reference 11 
or l4. In this case the lift-curve 
slope of the complete wing is as 
follows : 



m 2 - mi s x 2 
m 2 + mi E 


JT 

2 


(1 



+ 


_1_ 

m 2 



+ 


m 2 c 0 2 - 2c 0 ci(m 2 -m 1 ) 


m 2 (c 0 -c 2 ) 2 



(to) 


where S is the wing area and Ad, E, and K are the elliptic inte- 
grals (Appendix A) with modulus Jl- mi 2 and, in A 0 , the argument 
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Method of Computing Drag 

The drag on a flat -plate wing at an angle of attack is the differ- 
ence between the rearward component of the normal-force vector and the 
thrust due to leading-edge suction. Using linear theory to evaluate 
the first of these forces, we write 

D = a L - T (4l) 


The thrust T is calculated by integrating along the subsonic portion 
of the leading edge the streamwise component of the local suction force. 
If x,yi(x) are the coordinates of a point on the leading edge and m 
is the ratio of the inclination of the leading edge at that point to the 
tangent of the Mach angle, then, from reference 8, the contribution to 
the thrust at the point x,yi is given by 


dT 

dx 


m 


u ( x >y) P(yi-y) 


(42) 


where u Is the streamwise component of the perturbation velocity. 


Over the inside portion of the wings under consideration, m = mi, 

m i mi Uq x 

y x = — x, and u, near the leading edge, has the value — 

as In equation (12), from which is obtained 


dT Pit 2 r~-\ 2 ” /j.o\ 

— = — uq V 1 - mi x; x < c x (43) 

dx 2 

Then the total thrust over the inner portion of the leading edge is 

2 J H dx = | jru 0 2 ci 2 J~l- mi 2 (44) 

o 

If m 2 is greater than 1, equation (44) gives the total thrust. 
If m 2 Is less than 1, there is thrust along the outer portion of the 
wing as well. The velocity u In this region is found from equa- 
tion (20). In equation (20), the term containing G n remains finite 
as the leading edge is approached, so that it contributes nothing to 
the limit of the product in equation (42). The remaining term gives 
for the local thrust 


dT 

dx 



1 - m2 


XFrji(t) 


1 + m 2 


(45) 
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where 


P T (t) 


(l+mr) (t+mi) 


2E(k) 


mi 


K(k) 


12 


(W) 


k = 


/(1-mi) (t-mi) 
(l+mi) (t+mi) 


t = m 2 -(m 2 - mi) 


£1 

x 


While it is not feasible to integrate equation (45) analytically, it 
has been found possible to approximate the product xFm(t), within the 
accuracy of our computations, by a linear expression which, when inte- 
grated, gives 


-Pz 


dT 

dx 


2m 5 


dx 3? 


PUq* 


— 10.2 

l+m 2 


(c 2 - ci) 


Cl 


c 2 F t (To)- i (ci+c 2 )F T (m 2 ) 


(47) 


If a portion of the leading edge lies within the Mach cone from 
the bend in the opposite leading edge, as in sketch (e), there may be 
a considerable effect on the leading-edge thrust on that portion. 

Making use of the approximation of equation (23) to modify the strength 
of the leading-edge singularity in the pressure, retaining the linear 
approximation for xF T (t), and dropping the square of the small correc- 
tion contributed by equation ( 23 ), result in an additional term in the 
thrust: 


C 2 


x* 


dAT 

dx 


dx 


2m2 -> m 2 - mi / l-m 2 

— — PUo — p==== ( 

rtm i V 2(l+mi) \l+m 2 


3/2 


P T ( t o) 

L F T (m 2 ) 


c 2 +(c 2 -x*) 


yc 2 (c 2 -x*)F T (T 0 ) + 


*/Frp(m 2 ) 


f t( t q ) 

LF T (me) 


-,2 


c 2 -(c 2 -x*) 


cosh" 


2 Frp(Ti) 


c 2 F T (T x ) -(c 2 -x*)F T (m 2 ) 


m 


The drag is then calculated by substituting in equation (4l) the 
thrust given, by equation (44) and, if applicable, equations (47) and (48). 
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Numerical Results - Lift and Drag 


The lift and the drag due to lift were calculated for a large 
numher of wing plan forms , without , however, more than sampling the 
possible effects of the numerous parameters . The results are su mm arized 
in figures l4, 15 , and 1 6 . 

In figure llf, the ratio of inner to outer slope of the leading 
edge, the ratio of cx to the root chord, and the ratio of si to 
total semispan are fixed. Comparison is then made between the type B 
wing and a conventional sweptback wing of the same area, aspect ratio, 
and average sweep through a range of aspect ratios or, alternatively, 
through a range of Mach numbers, since any one curve can be interpreted 
as showing the effect of Mach number on a specific plan form. 


Sketch ( j) shows the method of deriving the conventional swept wing 
from the wing of type fe. The tip chord was kept unchanged for the com- 
parison because of the large loss in lift 
caused by the tip effect. The term 


"average sweep" is used in only a general 
sense, since no 
said to exist. 


exact equivalence can be 


Comparison Is also made with a 
type A wing having the same mx/m 2 
and Si/s ratios and approximately the 
same ci/c Q ratio. A single calcula- 
tion made with Ci/c 0 exactly equal to 
the value used for the type B wings 
showed this parameter to have little 
effect under the circumstances. 



Sketch (j) 


Curves showing the variation with aspect ratio of the lift and 
drag of rectangular and triangular wings are included for reference . 

For figure 1^, the aspect ratio 
was changed by the expansion of all 
the spanwise dimensions of a basic 
wing. In figure 15 the increase in 
aspect ratio is obtained by adding 
to a basic triangle triangular side- 
panels of increasing spanwise extent, 
starting at a fixed point back of 
the apex of the wing (see sketch (k)). 

Thus a family of cranked wings of 
type A is derived with fixed Ci/c Q 
and mi, but increasing m 2 , area. 



Sketch (k) 
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and aspect ratio. The resulting variations in lift and drag-rise factor 
are shown for two values of mi and two values of ci/c 0 . It should 
be noted that, since mi is kept fixed, the variation shown does not 
correspond to the effect of Mach number change. 

The curves for the sweptback wings in figures 1^ and 15 all show 
a more or less abrupt break when the leading edge coincides with the 
Mach lines, followed by a levelling off of the curves at higher aspect 
ratios. In the cases of discontinuously swept wings, the breaks appear 
when m 2 = 1, since mi was kept less than 1 throughout the calcula- 
tions. The effect of variations in the sweep angles of the leading 
edge is further shown in figure 1 6 , where the aspect ratio is constant, 
and only the shape of the leading edge is varied. The cranked wings 
considered in this figure are all derived from an aspect ratio 2 tri- 
angle by increasing the sweep of the leading edge at the root and adding 
compensating area at the tips (see sketch ( 2 )). The results are plotted 

against the ratio Ci/c D - the chordwise 
extent of the region of increased sweep - 
for selected values of mi, and cross - 
plotted for constant values of m 2 . 

DISCUSSION OF RESULTS 

The over-all characteristics of the 
cranked wings studied show no striking 
effects of the concentration of sweep 
near the root, at least as far as the 
lifting condition is concerned and omit- 
ting any consideration of the effect on 
the boundary-layer flow. The character- 
istics of the wings appear to be deter- 
mined primarily by the sweep of the 
outer portion of the leading edge, the 
increase in lift in these regions due to the more highly swept portions 
ahead being nullified by the decrease in lift in the center due to the 
small apex angle of the wing. In fact, the lift on the fully tapered 
(type A) wing is just about that on the triangular wing formed by the 
extension of the outer portions of the leading edge inward to the center 
line. Compared with triangular wings of the same aspect ratio, these 
wings show both lower values of the lift coefficient and higher values 
of the drag-rise factor. The calculations for the type B wings show a 
higher lift-curve slope than for conventional wings of the same average 
sweep and aspect ratio as long as the outer portions of the leading edge 
are subsonic, but also show higher drag due to lift. 

The use of a cranked leading edge has a marked effect on the distri- 
bution of lift. Comparison with the triangular wing having the same 
leading -edge sweep as that of the outer panels shows the pressure at 
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the rear to be essentially unchanged by the addition of the forward 
peak within the overlapped Mach cones from the bend in the leading edge., 
but to be reduced abruptly ahead of this region, approximately in pro- 
portion to the apex angle of the wing. The lift near the outer portion 
of the leading edge, on the other hand, is increased by the upwash from 
the forward area and, if the outer portion of the leading edge is super- 
sonic, the theory indicates infinite pressure along the outgoing Mach 
line from the bend. Somewhat the same effects may be anticipated behind 
the juncture of wing and fuselage in a complete airplane conf iguration. 

The center of pressure of the cranked wing lies behind the center 
of area if the tips are tapered to a point, but may lie ahead of the 
center of area if the tip chord is a considerable fraction of the root 
chord, because of the large loss of lift inside the tip Mach cone. In 
the range of Mach numbers investigated there was no movement of the 
aerodynamic center with change of speed. 

Any general conclusions regarding the value of the plan forms 
studied would have to take into account the possibilities offered in 
the way of a redistribution of thickness and the effect of the variation 
in sweep on the boundary -layer flow. The similarity of some of the plan 
forms to those of wing-body combinations suggests also that comparisons 
should really be made with such complete configurations, taking into 
consideration questions of total drag, weight and useful volume. These 
questions, though of prime importance, are beyond the scope of this 
report. 

Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., March l8, 1953- 
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APPENDIX A 
SYMBOLS 
General 

free -stream velocity- 
free -stream Mach number 
density of air 
angle of attack, radians 
a/ M 2 - 1 

horizontal perturbation velocity, positive downstream 

vertical perturbation velocity, positive upward 

streamwise component of the perturbation velocity along 
the center line of a flat triangle at an angle of attack, 
with semiapex angle less than Mach angle (eq. (3)) 

upwash of a flat triangle at an angle of attack a, with 
semiapex angle less than Mach angle (eq. (l)) 

lift 

drag 

thrust 

Ifu 

coefficient of lifting pressure, — 

coefficient of lifting pressure on flat triangle at angle 
of attack a, with semiapex angle less than Mach angle 
(eq. (12)) 

coefficient of lifting pressure added by presence of each 
outer panel 

coefficient of lifting pressure added by presence of tip 
Wing Dimensions 

wing area 


wing semi span 
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root chord (fig. 2) 

streamwise distance from apex of wing to bend in leading 
edge (fig. 2) 

streamwise distance from apex of wing to tip of leading 
edge (fig. 2) 

spanwise distance from wing center line to bend in leading 
edge 

P times the slope of the inner segment of the leading edge 

P times the slope of the outer segment of the leading edge 

Coordinates 

x,y Cartesian coordinates, positive downstream and to the 

right of the flight direction, respectively 

Xi,yi coordinates of point on the leading edge (eqs. (8) and ( 9 )) 

s coordinates of intersection of ray from Xi,yi with tip 

chord 

y*(x) y expressed as a function of x for points along the 

inboard Mach line from the point c 1; si 

x*,y* coordinates of intersection of inboard Mach line from ci,si 

with the leading edge of the left-hand panel 

t Q P times the slope of a line from the apex of the wing 

ti P times the slope of a line from xi,yi to the point at 

which the pressure is being calculated (eq_. (4)) 

t t Q of the point at which the pressure is being calculated 

r 0 ,r lf t 2 ,t 3 limiting values of t Q (t q given by eq. (10), t 1 by 
eq. (25), others defined where used) 

t* P times the slope of a line from x*,y* 

t m P times the slope of a line from ci,Si (eq. (7)) 

t s P times the slope of a line from c 2 ,s (eq. (34)) 


c o 

Cl 

c 2 

si 

m x 

nig 
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T q 0 times the slope of the line joining the apex and the 

leading-edge tip (eq. (39)) 

Ti,T 2 ,T 3 0 times the slope of the lines from x x ,yi to (c 2 ,-s), 
(c Q ,-s), and (c Q ,s), respectively 


Mathematical Symbols 


P 


P1.P2 
Qi*Q2 
Gn; S 

TL 


function of t x proportional to lifting pressure in 
upwash-canceling solution (eqs. (l4) and ( 17 )) 

functions entering into P for m 2 <l (eqs. ( 15 ) and (l6)) 

functions entering into P for m 2 > 1 (eqs. (l8) and ( 19 )) 

functions entering into Z£p for m 2 <l (eq. (20)) 

function of t Q proportional to the total lift induced by 
an upwash-canceling solution 

function of t entering into the calculation of the 
leading -edge thrust (eq. (46)) 


k 


modulus of elliptic integrals, defined where used 


k' 


complementary modulus ( J 1 - k 2 ) 


<p,0, or ijr argument of elliptic integrals, defined where used 

F(9»k) incomplete elliptic integral of the first kind of argu- 

ment cp and modulus k 


E(cp,k) 


incomplete elliptic integral of the second kind of argu- 
ment cp and modulus k 


K,K(k) complete elliptic integral of the first kind; that is, 

K = F (f- k ) 

E,E(k) complete elliptic integral of the second kind; that is. 



Jacobian zeta function; Z(cp,k) =* E(cp,k) 
(KZ(cfk) is tabulated in ref. 11.) 


. IM. F (cp,k) 
K(k) 


Z 
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A 0 function used in evaluating elliptic integral of third 

kind, circular case; 

Ao(cp,k) = ||K(k)E (9,k') - [K(k)-E(k)]F(cp,k')j 
(tabulated in refs. 11 and l4) 
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APPENDIX B 
FORMULAS FOR LIFT 

For m 2 < 1, 

IL(to) = -2s 2 Pi(t 0 ,Ti) + , ^i -( s +o g t 0 ) 2 %(t Q ,Ti) - 

m 2 +t 0 


(c 0 -Xi ) 2 (m 2 -T 2 ) 2 Pi(to,T 2 ) -(t 0 -T 2 ) 2 P 2 (t 0 ,T 2 ) ~ 



The terms in Ti and T 2 vanish when Ti and T 2 equal "1* 
cone from xi,yi does not include the far tip, Ti and T 2 are less 
than -1 and the inapplicable terms are imaginary. For wings of type A, 
with zero tip chord, Ti equals T 2 and T 3 equals m 2 , and the formula is 
accordingly simplified. 

For m 2 = 1, 






;-t 0 / 2 (mg-tp) 

2 +t 0 V m 2 ( 1 - 212 ) (l-to) 


sihh -1 


gm 2 (l+T x ) 

(l+m 2 ) |m 2 +Ti| 



NACA TO 2959 


31 


For 

IL (^0) 


For 

3 L(t 0 ) = 


m 2 > 1 , 


= (c 0 -xi ) 2 


(m 2 -T 2 ) 2 Qi(t 0 ,T 2 ) - (t Q -T 2 ) 2 Q 2 (t 0 ,T 2 ) - 


(m 2 -to) f sin -1 _ sin-i /t 7) .. . cos -1 + 

V V 2 ' A/ m 2 (m 2 -l) J m 2 -T 3 

(Vfa)f 1or >0(1^3) + /T 3 (l-t Q ) - 
Vt 0 (l-t 0 ) 7 it 0 -T 3 | 


m 2 = 00 } 


(c 0 -ci ) 2 


(t 0 -2T 2 )sin -1 



(t 0 -T 2 ) 2 Q 2 ("tQ^T 2 ) + 


\(j T 3 (l-T 3 )- /I^?) + ( 2 T 3 -to- 0 sin-1 t 3 + 

(tp-Tg ) 2 iog /t Q (l-T 3 ) +y T 3 (l-t 0 ) - 

y^> j [t 0 -t 3 | 
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Type A 



Figure 1.— Plan forms typical of the two classes considered. 
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Figure 2.— Cartesian coordinates and plan— form parameters. 
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Figure 3.— Upwash field for flat triangular wing at an angle of attack and 
.in moderately supersonic flight (Mach lines ahead of leading edge). 
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Figure 7.— Plan form of wing for which lift distributions were calculated, 

and Mach line pattern at M = 1.4l4. 
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Figure 8.— Spanwise lift distribution at various streanwise stations of wing shown in figure J, 

with M = 1.4l4 (m 2 <l). 
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Figure 11.— Variation of aerodynamic— center location with. Mach number 
for wing of figure 7, and for type A wing with same leading edge. 
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Figure 13 . ~ Method of correcting lift of elementary upvash-cancellation 

field for tip effect. 
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Wing $ c f/ Co * 0.294 
Wing A, c t/c. = .294 



%M£ A r*. f 

ntng A, x 'i/c* x ovo 


Conventional wings, 
same average sweep 
as B (see text) 


/9x Aspect ratio 

(a) Lifixiurve eloper 

Figure Vfijri&"bion of lift curve slops 0 anrl Ixsig 3riss factor Qr\^/ r with aspect ratio 

and Mach number. (Aspect-ratio increase obtained by increasing all spanvise dimensions.) 
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with m 1 ^m 2 and chordwise 
X aspect ratio = 2. 





